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I. INTRODUCTION 



The advent of quantum algebras precipitated many new results in the area of in- 
tegrable models. Likewise their supersymmetric counterparts quantum superalgebras 0-|6j 
facilitate systematic treatments of integrable models which accommodate both bosonic and 
fermionic degrees of freedom. An important subclass are those which may be interpreted 
as describing systems of correlated electrons for their obvious physical applications in con- 
densed matter physics. Among these are the supersymmetric (SUSY) t — J model and 
supersymmetric generalizations of the Hubbard model |],[TTJ • Particularly, it would be ben- 
eficial if these models could provide some knowledge of the transitions between the metallic, 
insulating and superconducting phases. It has long been known that an insight into these 
properties may be gained by studying the effects of the boundary conditions on such models 
(e.g. see jIT| , [l2| ). Subsequently several authors have studied electronic models with twisted 



boundary conditions ||iq-|lT|. 

Integrable models with twisted boundary conditions may be formulated within the frame- 
work of quantum (super) algebras. This notion of twisting is more general than that usually 
treated in electronic models whereby a twisted boundary condition is thought of as the in- 
troduction of a phase factor into the periodicity of the model. Here our twisted boundary 
conditions correspond to more general transformations of the local states and have their ori- 
gins in the underlying symmetry of the model [ 20| . For a given i?-matrix R(u) any matrix 
M satisfying [M1M2, R{u)] = allows an integrable model to be constructed with the trace 
over the auxiliary space weighted by M. This is what we will refer to as a twisted boundary 
condition. 



The work of Reshetikhin |21] on multiparametric quantum (super) algebras permits a 
natural way to construct integrable models also dependent upon additional free parameters. 



In fact the model of Perk and Schultz p2| may be formulated within this framework ]23 



We will make it apparent that there exists under suitable conditions a integrability pre- 
serving mapping between these multiparametric models and those with twisted boundary 
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conditions. Some particular cases have been studied previously |24j-g7fl. We show that these 
methods may be applied to any model with an underlying quantum (super) algebra symme- 
try. By establishing this correspondence between models with twisted boundary conditions 
and multiparametric models, it is reasonable to expect that the multiparametric solutions 
may provide suitable test models for describing the various phases associated with correlated 
electron systems. 

As examples we will consider the cases of the SUSY t— J model [§] and the SUSY U model 
ig. As well as giving the Hamiltonians for the integrable multiparametric generalizations 



of these two models, we have also determined the corresponding Bethe ansatz equations 
which provide the starting point for an investigation into their thermodynamics. 

The paper is organized as follows. In section II we present the general construction of 
multiparametric spin chains and their relation to models with twisted boundary conditions. 
In section III we illustrate how our general formalism applies to construct multiparametric 
versions of the SUSY t-J and U models. The Bethe ansatz equations of the models are also 
obtained. A summary of our main results is presented in section IV. 



II. GENERAL CONSTRUCTION 

Using Reshetikhin's construction [plj for multiparametric quantum algebras, it is 
straightforward to obtain the associated multiparametric quantum spin chain. Here we 
demonstrate that under appropriate constraints these models may be transformed to quan- 
tum spin chains with twisted boundary conditions; i.e. the additional parameters arising 
from Reshetikhin's construction may be mapped to the boundary. 

Let (A, A, R) denote a quasitriangular Hopf (super) algebra where A and R denote the 
co-product and _R-matrix respectively. Suppose that there exists an element F G A® A such 
that 

(A <g> I)(F) (I® A)(F)= F13F12, 

F12F13F23 = F23F13F12, F12F21 = I (1) 
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Theorem 1 of [21] states that (A, A , R ) is also a quasitriangular Hopf algebra with 



co-product and i?-matrix respectively given by 

A F = F 12 AF 21 , R F = F 21 RF 21 . (2) 



In the case that (A, A, R) is an affine quantum (super) algebra we have from [21| that F 
can be chosen to be 

F = exp ^2 ( H i ® Hj - Hj ® H) fa (3) 

i<j 

where {Hi} is a basis for the Cartan subalgebra of the affine quantum (super) algebra and 
the fa, i < j are arbitrary complex parameters. For our purposes we will extend the Cartan 
subalgebra by an additional central extension (not the usual central charge) H Q which will 
act as a scalar multiple of the identity operator in any representation. 

Suppose that 7r is a loop representation of the affine quantum superalgebra. We let R(u), 
R F (u) be the (super)matrix representatives of R and R F respectively, which both satisfy 
the Yang-Baxter equation 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 {v)R n {u)Ri 2 {u - v). 

For the supersymmetric case it is necessary to impose the multiplication rule 

(a®b)(c®d) = (-l) [b][c] (ac®bd) 

for homogeneous supermatrices a, b, c, d where [a] = if a is even and [a] = 1 if a is odd 
p8| . However with an appropriate redefinition of its matrix elements R(u) satisfies the 



usual (non-graded) Yang-Baxter equation (e.g. see |p9| ). Throughout we adopt this latter 
convention. 

If R{u)\ u=Q = P with P the permutation operator then R F (u 
(pi). We may construct the transfer matrix 



result of 

u=0 



t F (u) = str (vr^ +1 ) (l®A F ) <) 

= str (R F N (u)R F {N _ 1) (u)....R F 1 (u)) (4) 
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where A^ is defined recursively through 

A F = (J® J....® A F ) A^_ x 

= (A F ®/....®/) A^_!. (5) 

The subscripts and 1,2,...,N denote the auxiliary and quantum spaces respectively and 
stro is the supertrace over the zeroth space. From the Yang-Baxter equation it follows 
that the multiparametric transfer matrices t F (u) form a commuting family. The associated 
multiparametric spin chain Hamiltonian is given by 

^=( t »)- l A t »| 

u=0 

AT-1 

= E ^+1 + hm (6) 



with 



h F = 4- PR F {u) 

du 



u=0 



Through use of ([!]) we may alternatively write 

t F (u) = str (tt^ 1 ) (/ ® J^v) [(/ ® A Ar )(F 10 i? 01 F 10 )] (J ® J^ 1 ) 



with 



Jat — Gn~iGn-2----Gi, 

G% — F iN Fi( N _iy...Fi(i +1 y (7) 



We now define a new transfer matrix 



t(u) = J^H F {u)J N 



= str (vr®^ 1 ) (/ ® A N ) (F W R 01 F 10 )) (8) 

where we have employed the convention to let F denote both the algebraic object and its 
(super)matrix representative. Through further use of ([!]) we may show that 



t(u) = str (F 10 F20....F N0 R 0N (u)R 0iN _ l) (u)....R 0l (u)F 10 ....F 1 
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and the associated Hamiltonian is given by 
H = (tin))- 1 ±t(u) 



u=0 

iV - 1 > 2 



^ h*i t i+i + (Fn(N-1)--~Fni) hjsfi (F 1A r....F( i v-i)Ar) (9) 



i=l 



where 



h = ^-PRiu) 
du 



u=0 

The above Hamiltonian describes a closed system where instead of the usual periodic 
boundary conditions we now have a more general type of boundary condition. The boundary 
term in the above Hamiltonian is a global operator; i.e. it acts non-trivially on all sites. 
However we can in fact interpret this term as a local operator which couples only the 
sites labelled 1 and N. It can be shown that the boundary term commutes with the local 
observables h^+i for i ^ 1, N — 1. This situation is analogous to the closed quantum 
(super) algebra invariant chains discussed in p0| . 

From the above construction we may also yield models with twisted boundary conditions 
by an appropriate choice of F. Recall that we extend the Cartan subalgebra by the central 
element H . Let this element act as cl in the representation it where c is some complex 
number. If we now choose 0^ = for % ^ in the expression ([3D the matrix F factorizes as 
F = M± 1 M 2 with 

M = exp c^oiHi^j 

and / is the rank of the underlying quantum (super) algebra U q (g). Using the fact that the 
i?-matrix satisfies 

[R(u), I ® Hi + Hi ® I) = 0, i = 1, 2, I 

tells us that 

[R{u), M X M 2 \ = 0. 
In this case the Hamiltonian (|9]) reduces to 

AT-l 

H=J2hi+i + M? N h m Mr 2N (10) 
1=1 

which is precisely the form for a system with twisted boundary conditions (see |20[ ) . 



III. EXAMPLES 



In this section we illustrate how our formalism applies to construct a multiparametric 
version of the SUSY t-J model and the SUSY U model fTU |. Both models are gl{2/l) 
invariant and their formulation through the quantum inverse scattering method can be 
found in || and |JT| , respectively. The first model describes electrons with nearest-neighbor 
hopping and spin exchange interaction on a chain, while the second can be considered an 
extension of the Hubbard model with additional pair-hopping and bond-charge interaction 
terms. These models are of interest because of their possible connection with high-Tc su- 
perconductivity. In order to turn our discussion more general, we will in fact handle with 
their anisotropic or q-deformed versions [ 32| , |3^j3^]. Of course, in the rational limit q — > 1 
all results reduce to their corresponding isotropic ones. 



A. The supersymmetric t-J model 



We begin by introducing the multiparametric U q (g£(2/1)) i?-matrix, which in terms of a 
generic spectral parameter x and a deformation parameter q reads 



where 



a = xq — 



R F (x) 



xq 
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(12) 



and ti,t 2 ,t3 are independent parameters written in terms of 0oi,0o2,0i2 as 

t x = exp(-20 O i + O2 + 12 ) 

t 2 = exp(-0 O i + 002 + 012) (13) 
t 3 = exp(0 ol - i2 ) 



The above matrix was already presented by Perk and Schultz [£4] when studying a 
multicomponent generalization of the six-vertex model. 

Next we construct the transfer matrix t F (x) according to eq.([|) from which we find the 
associated multiparametric Hamiltonian (see eq. (|) on a one-dimensional periodic lattice 

JV-l 

i=i 



where 



— ~~ (T2 C iT C *+lT + ^2 c i+lT C *T + T2 C i| Ci +U + ^3 C !+U C i|) 



6 2 6 3 



- 2 



+ t?5r5i ! + 2 cos 7 (5/5^ - 



+ isin(7)(nj - n i+ i) - z sin(7)(n i S' J 2 +1 - S?n i+1 ) - cos 7^ + 2 cos 7 (14) 

Above c£2's are spin up or down annihilation (creation) operators, the S^'s spin matrices, 
the Hi's occupation numbers of electrons at lattice site i and 7 is the anisotropy parameter 
(q = e %1 ). A similar version of a multiparametric SUSY t-J model has already been discussed 
in ref. [fjq] . Notice that here it emerges systematically from our general construction. By 



setting ti,t 2 ,t s — > 1 in eq. (|14]) , the usual terms of the anisotropic SUSY t-J model |32] can 
be recovered. 

The Hamiltonian (|14"|) can be exactly solved through the algebraic nested Bethe ansatz 
method. This procedure is carried out in two steps and the Bethe ansatz equations are given 
by 



+ 2(N-M 2 ) + 2M2+-2M2 ( a ( X k > ) \ TT a ( x P / x< k^ ) H x ^/ X i^) TT ^ I X k^) _ , _ -> M 

1 2 3 1 77 (i)J 'Li (i)/ (iK 7 Wi (i)7 11 7 (2)/ (1)7 ~ ' 1, . . . iWi, 

V K4 ) y *=i K4 74 ) a(4 74 )i=ia(4 74 ) 

, M X / (2) / M 2 / (2) / (2K , / (2) / (2)x 

-2(JV-M 1 ),2(JV-M 1 ),2Mi/ i \AT a TT / X ' J TT /Xfc I A k 13. > - 1 fc - 1 M (\ <\ 

k h h M) S M K*?»/4 2 ') Af) " " " ' ■ ■ ■ 2 ( ) 



where x 



M, 



1,2; k = l, 



M m ) denote the Bethe ansatz parameters, N is the number 



of lattice sites, Mi is the number of holes plus down spins and M2 is the number of holes. 
We see from (ITS!) that the additional parameters £1, £2, £3 have the meaning of external fields 
(see e.g. 



36,3 



Following the approach presented in the previous section, we perform the transformation 
(§) and then set t 3 = ^ (or 0i 2 = , see eq. (|I3D) in order to find the anisotropic SUSY t-J 
model with twisted boundary conditions ( JTH ) 



2V-1 



Nl, 



8=1 



where = lim{ tlit2)t3 _»i} and 



27V I | 
l 2 



tl 



2N 



c ni c h + 



4 



2/V 



+ i sin(7)(riiv — ^1) — i sin(7) (njyS'i — S z N n\) — cos 7?^ + 2 00S7 



(16) 



B. The supersymmetric U model 

Let us now construct a multiparametric version of the anisotropic SUSY U model, which 



has been proposed recently as a new integrable model for correlated electrons (see ref. fl33|j34 
for more details). 

We begin by recalling the trigonometric R-matrix associated with the one parameter 
family of four-dimensional representations of U q (g£(2/1)) 

R(x) = PR(x), 

„x „2a 1 „x+2a+2 

R( x ) = ? — q —Pi + P 2 + — —r-ur p z- (17) 

Here x and q are, respectively, the spectral and deformation parameters and a is a free 
parameter which arises from the underlying representation. P is the permutation operator 
and Pi, z = 1,2,3 are projectors whose explicit form can be found in fl34 |. 
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We find the corresponding multiparametric i?-matrix 



R F (x) 
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(18) 



where t 1 , t 2 and t 3 are independent parameters also given by eq. ( |T3D and "*" denote the 
elements of the R- matrix (|TTD, which can be obtained from the projectors given in [M\. We 
do not write them explicitly here since we will not need them later. Notice that here, in 
contrast to the previous case (see eq. ([TT|)) , the new parameters £1,^2, £3 occupy also non- 
diagonal positions. This is a peculiarity of higher representations and can also be verified 
for other higher spin models ( e. g. spin 1 XXZ chain). In fact, our prescription for the 
element F of the multiparametric R matrix is particularly interesting in these cases, where 
it is not obvious how to construct R F . 

Next we construct the transfer matrix t F (x) according to eq.(||) from which we find the 
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multiparametric version of the anisotropic SUSY U model (see eq.( 

N-l 

H F =T,hf, +1 + h F N1 , 
i=i 



where 



Ki+i = -^ + iA-ri- l v + TH-ri- l v-} rH+n ~ h.c. 
-p4c i+ uh^^] niT [-^^ + ] ni+1T - h.c. 

+ { a \ q 1 £p c ii4] c i+^ c i+n + h.c. 

+<f +1 (n iT + n n - 1) + cr^+iT + n i+u - 1) (19) 



and £ = ^§ , 77 = ^ , p = 72W , = sgnicAq^ ( ^rr^ )'- In the h.c. terms one should 
notice that the parameters ti — > t' 1 ^ = 1,2,3. Equation ( 19 ) is a generalization of the 



anisotropic SUSY U model |34|, which can be recovered in the limit ti,t 2 ,t 3 — > 1 . 

This model can be exactly solved by means of the algebraic nested Bethe ansatz method 
and the Bethe ansatz equations are given by 



W^a+l „-l\JV ^2 ™W _ ~( 2 )„2 



k = 1, . . . JVi 



,2(N-N 2 ).2N 2i .-2N 2 ( _TT 

1 2 3 ^ 4V«-i J "M,(4"-xf, 

, JVi _(2) J V 2 _ (2) , (2) 2 

*i *2 *3 11— 711 f^-ll - (2I ftTT ' fc-l,...iV 2 (z(Jj 



q(x\ -xl>) x k > - x) >q 2 

where x k (m = 1,2 ; k = l,...,N m ) are the Bethe ansatz parameters, N\ is the total 
number of spins and N 2 is the number of spins down. 

According to the approach presented in the previous section, we perform the transfor- 
mation (||) and then choose t 3 = |^ in order to find the anisotropic SUSY U model with 
twisted boundary conditions, which yields 

N-l 

H = } J h i)i+ i + h N1 , 



where /i» i+i denotes the local terms of the anisotropic SUSY U model [34| and 
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h NA = -tf c+ T c 1T (-z/ + )^(-^) nu " h.c. 
-tl N c+ l Ci l (-v-) nN i{-v + ) n n - h.c. 
+ [a]g 1 (t 1 t 2 )~ 2N 4- 4 4" T ci T c u + h.c. 
+ [a]q 1 ( n mn Nl +n iT n u ) 

+q a+1 (n m + n Nl - 1) + g-^niT + n n - 1) (21) 



IV. CONCLUSIONS 

In this paper we have demonstrated a correspondence between multiparametric spin 
chains and models with twisted boundary conditions in the expectation that this connec- 
tion will provide further insight into the description of phase transitions of such integrable 
systems. Our approach can be applied to any model with an underlying quantum (su- 
per) algebra symmetry. We are particularly interested in models which describe systems of 
correlated electrons and have studied the SUSY t — J and U models as examples. 

Another important class of integrable models are those associated with the Temperley- 
Lieb algebra. In Zhang proposes a systematic method to generate multiparametric 
extensions of these models. It is possible to adapt the techniques employed in this paper to 
establish a mapping from models based on the Temperley-Lieb algebra with twisted bound- 
ary conditions to associated multiparametric generalizations. With respect to correlated 
electron systems, an example based on the Temperley-Lieb algebra has been described in 
39| , |40[| to which this procedure can be applied. 
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